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ABSTRACT 

A magnetic field is force-free if there is no interaction between the magnetic 
field and the plasma in the surrounding atmosphere i.e., electric currents are 
aligned with the magnetic field, giving rise to zero Lorentz force. The com- 
putation of various magnetic parameters such as magnetic energy (using virial 
theorem) , gradient of twist of sunspot magnetic fields (computed from the force- 
free parameter a) and any kind of extrapolations, heavily hinge on the force-free 
approximation of the photospheric sunspot magnetic fields. Thus it is of vital im- 
portance to inspect the force-free behaviour of the sunspot magnetic fields. The 
force-freeness of the sunspot magnetic fields has been examined earlier by some 
researchers ending with incoherent results. The accurate photospheric vector field 
measurements with high spatial resolution are required to inspect the force-free 
nature of sunspots. For this purpose we use several vector magnetograms of high 
spatial resolution obtained from the Solar Optical Telescope/Spectro-Polarimeter 
(SOT/SP) aboard the Hinode. Both the necessary and the sufficient conditions 
for force-freeness are examined by checking global as well as local nature of equi- 
librium magnetic forces over sunspots. We find that the sunspot magnetic fields 
are not much away from the force-free configuration, although they are not com- 
pletely force-free on the photosphere. The umbral and inner penumbral fields 
are more force-free than the middle and the outer penumbral fields. During their 
evolution, sunspot magnetic fields are found to maintain their proximity to force- 
free field behaviour. Although a dependence of net Lorentz force components is 
seen on the evolutionary stages of the sunspots, we do not find a systematic re- 
lationship between the nature of sunspot magnetic fields and the associated flare 
activity. Further, we examine whether the fields at the photosphere follow lin- 
ear or non-linear force free conditions. After examining this in various complex 
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and simple sunspots, we conclude that, in either case, the photospheric sunspot 
magnetic fields are closer to satisfy the non linear force-free field approximation. 

Subject headings: Sun: photosphere. Sun: sunspots, Sun: surface magnetism 



-4- 



Introduction 



fll954h 



'he concept of force-free magnetic fields first came into existence wlien iLiist fc Scfiliiter 
pointed out that the cosmic magnetic fields often follow a condition given as, 



Equation 1 represents zero 



1979; 


Low 


1982 


)• 



(V X B) X B = 0. 



orentz force and the field following this equation is 



Known 



Chandrasek. 



Chandrasekharl ( 1956a 



lar fc Kendal 



1957 



Chandrasekhar 



1961 



Parker 



b|) explained many characteristics of the force-free 
field equation including its solutions in axisymmetric as well as in non-axially symmetric 
cases for a = constant. Man y of the c harac t eristics of force- free behaviour in the so lar 



atmosphere were discussed by 



1973 



1974 . 



1977 



1978 



19801 ) 



Woltierl ( 



19581 



Mo: 



Low &: Nakagawal (119751 ): 



odenskii ( 



1969) 



Molodenskvl ( 



Lowfc Flvej (120071 ): 



Parked ( 11989 



Low 



19901 ) etc. 

We can rewrite Equation 1 as, 



V X B 



aB 



(2) 



which denotes that the electric current density J (= l/47r V x B) is proportional to the 
magnetic field B. Here a is a constant and some function of position. If we take divergence 
of the Equation |2l as a result of Maxwell's equation 



V-B = 



we get the following 



B ■ Va 



0. 



(3) 



(4) 



Equation 4 shows that the force-free parameter a remains constant along any field line. 
However, the a can vary across the field lines. The z component of Equation [2] allows us to 
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compute the distribution of a on the photosphere (z=0): 



a 



,dBy OB, 



dx dy 

Three cases may arise: (i) a = everywhere, i.e., there is no electric c u rrent i n the 



(5) 



atmos p here resu l ting in a po t ential field (for d e tails 



Sakurail (11982 



19891); 



Regnier fc Priest 



1^ 



1989 



i near force-free field state (see e.g 



please see, 



Schmidt 



(11964f l: 



Semel 



20071)1 (ii) a = constant everywhere, 



van Ballegooijen fc Cranmer 



2010 



Nakagawa &: Raadu 



1972 



Gary 



1989 



Sakurai 



etc) which is not always yalid a n d (iii) a ya ries 



spatially, i.e . , nonlinear 



Amari et al. 



1997, 



1999 



'orce-free magnetic field (see e.g.. 



2006 



Wiegelmann 



2004: 



Sakurai 



Regnier fc Priest 



1979 



2007 



Low 



1982 



Schriiver et al. 



2006 . 



2009 



2008 



De Rosa et al 



2009 



Wheatland fc Regniei 



Metcalf et al 



3 



2008 



Mackay fc van Ballegooijen 



2006b 



20091 . etc). This is the most common state expected (iLow 



19851 ). However, the vector magnetograms with high spatial resolution are required to 



confirm the non-linear force-free state of sunspot magnetic fields. 

The potential field approximation of sunspot fields does not appear valid due to th e 



presence of 



Tiwari et al. 



■ ocal ve rtical current distributed over the sunspots (see, e.g. 



Su et al. 



20091 : 



2009bl . and references therein). The linear force- free field approx i matio n has 



Moon et al. 



(12002) using 



been found closely valid in a sunspot NOAA AR 5747 studied by 
relatively lesser resolution data. In earlier works, perhaps the insufficient resolution of 
the data obscured the conclusions about the validity of linear/non-linear force-free field 
approximations. In the present work, we also inspect the validity of linear or nonlinear 
assumptions after examining the force-free nature of sunspot magnetic fields by using 
high spatial resolution photospheric vector magnetograms obtained from Solar Optical 
Telescope/Spectro-Polarimeter onboard the Hinode. 

The general belief is that the photospheric field is not force-free due to high plasma 
P (ratio of plasma pressure to magnetic pressure) and is force-free in the corona due to 
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low plasma /3 (see e.g., iGaryl 1200 ll ). However this has become possible to verify on the 
photosphere with the advent of the vector magnetograms with high spatial resolution. 
Observations of vector magnetograms in the chromosphere will also make it feasible. 
However, such data are yet to obtain. Spatial variations of plas ma and magn e tic pr essures 



Solanki et al. 



(119931) and 



over a sunspot obser ved in infrared lines have been studied by 
Mathew et al.i (|2004) . The fields of the umbra were found more force- free than that of the 
penumbra. Here, we make an attempt to find the proximity of the photospheric field to 
force-free n ature by using the hi gh re solution vector mag netograms obtained from Hinode 



(SOT/SP). 



Metcalf et al 



(jl995|) and 



Moon et al. 



(120021 ) have already performed some 



analysis in this regard using ground based observations. By using the obse rvations of 



Stokes parameters at six wavelengths within the spectral line of Na I A5896, 



Metcalf et al. 



(119951 ) concluded that the magnetic field for an active region (AR) NOAA 7216 is not 



force-free in the photosphere whereas it becomes force-free a p proxi mately 400 km above 



Moon et al. 



(I2OO2I ). by examining the 



the photosphere in the chromosphere. Whereas 
observations of Fe I doublet A6301.5 and A6302.5 for 12 vector magnetograms of three 
active regions, found that the photospheric magnetic fields are not so far from the force-free 
na ture as conventionally regarde d. Another rigorous study in this regard was performed 



by 



data. 



Georgoulis fc LaBontd (12004f). W ith some limitations of their techniques and quality of 



Georgoulis &: LaBontd (l2004l ) concluded that the photospheric active region magnetic 



fields are not force- free. 

All the above results are not well in agreement. One of the reasons may be the 
inadequate resolution of their data sets. The accurate photospheric vector magnetic field 
measurements with high spatial resolution are needed to fully examine (both necessary and 
sufficient conditions) the force-free nature of sunspots. We have used 19 high resolution 
vector magnetograms obtained from the Solar Optical Telescope/Spectro-Polarimeter 
(SOT/SP) aboard the Hinode to inspect the force- free behaviour of the photospheric 
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sunspot magnetic fields. The effect of polarimetric noise present in the data o 



Dtained from 



SQT/SP does not affect m uch in derivation of the magnetic field parameters (ITiwari et al. 



2009a : 



Gosain et al. 



20101 ) 



We have further analyzed total 60 vector magnetograms of four active regions: two 
complex and fiare-productive, namely NOAA ARs 10930 and 10960, and two simple and 
relatively quiet, namely NOAA ARs 10933 and 10961. Total of 18 vector magnetograms 
for the AR 10930, 12 for the AR 10933, 15 for the ARs 10960 and 10961 each, have been 
analyzed to look at the variations in the force-free behavior during their temporal evolution 
on the one hand and to find any relationship with the flare activity on the other hand. 
We also investigate whether linear or non-linear force-free field approximation is closer to 
validity in the sunspot magnetic fields as observed from the SOT/SP. 

The paper has been organized in the following way. The methods to test the local and 
global force-freeness of the photospheric sunspot magnetic fields has been illustrated in the 
following section (Section 2). In Section 3, we describe the data sets used and methods of 
its analysis. Section 4 portrays the details of the results obtained. Finally in Section 5, we 
present our conclusions with appropriate discussion. 



2. Conditions for the Force-freeness of the Magnetic Fields 

Before considering the parameters which depend on the force-free parameter a and 
perform the modeling that are based on the force-free assumption, as discussed earlier, 
it is of vital importance to examine whether the photospheric magnetic fields are really 
force-free, and such approximations are realistic. For this purpose we need to examine 
both global and local nature of Lorentz forces, i.e., both the necessary and the sufficient 
conditions. 
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2.1. Necessary condition 



We use the fo 



Molodensky 



1974 : 



lowing well known prop erty of the Lorentz force (IMolodenskij Il969 



Parker 



1979 



Low 



19851 ) that the Lorentz force can be written as the 



divergence of the Maxwell stress tensor: 



in 



5ij + B^Bj/An 



(6) 



We can convert the tot al force over a volume in to the total Maxwell stress exerted upon the 



boundary surface (e.g., 



Parker 



1979 



Low 



19851 ). Under the assumption that the magnetic 



field above the plane z=0 (photosphere) falls off enough as z goes to infinity, the net L orent z 



force in the hal 



1984 



Low 



-space z>0 is just the Maxwell stress integrated over the plane z=0 (lAly 



19851 ). Thus the components of the net Lorentz force at the plane z=0 can be 



expressed by the surface integrals as follows: 



— / ByB^dxdy 



^ j {Bl -Bl- Bl)dxdy 



(7) 



(9) 



where F^., Fy and F^ are the components of the net Lorentz force. According to 
the necessary conditions for any magnetic field to be force-free are that 

I F^ |< Fp, I Fy |< Fp, and \ F, |< Fp. 



Lowl fll985f ) 



(10) 



Where Fp is the force due to the distribution of the magnetic pressure on z=0, as given by, 

= ^ / (Bl + Bl + Bl)dxdy. (11) 

Thus the inequalities given by Equation [10] provide a test for whether a force-free 
field exists at the boundary plane z = (photosphere). The criterion for the fields to 
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be force- free is that the Equations 7 , 8 and 9 normahzed with respect to Fp s 



be much less than unity (Lo 



3wl 



lould each 



19851 ). It was discussed by 



Metcalf et al 



fll995[ ) that the 



magnetic field is completely force-free if the aforementioned ratios are less or equal to 0.1. 



It is to be noted, however, that the inequalities given by Equation [TO] are only necessary 
condition for the fields to be force-free. The reason for this is that some information is lost 
in the surface integration in Equations 7, 8 and 9. 



2.2. Sufficient condition 



We can infer the sufficient condition for the force-freenes s of sunspot magnetic fields 



from the distributions of vertical tension force over sunspots (jVenkatakrishnan &: Tiwari 
20101 ). In a force- free case the tension force will balance the gradient of magnetic pressure 
demanding for zero Lorentz force. The usefulness of the te nsion force has not foun d 



much attention earlier in the literature but for few studies ( Venkat akrishnan 



Venkatakrishnan et al.l Il993 



Venkat akrishnan fc Tiwari 



I Venkat akrishnan fc Tiwarii |201Q : iTiwari 



1990a 



2010|). Recently 



(I2OIOI ) pointed out the utility of tension force as a diagnostic of 



dynamical equilibrium of sunspots. 



Let us begin with the equation of magneto-hydrostatic equilibrium 



(V X B) X B/47r - Vp + pfif = 



(12) 



where first term is the Lorentz force, second term is the gradient of gas pressure and the 
last term is the gravitational force onto the plasma. Lorentz force (let us call it as F) can 
be split into two terms as, 

„ (B.V)B V(B.B) 



The first term in the right hand side in the above equation is the tension force (let us call it 
as T). The second term represents the gradient of the magnetic pressure i.e., the force due 
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to magnetic pressure. We can expand the vertical component of the tension force in terms 
of the derivatives of the transverse components of the magnetic fields as follows: 



An ox oy 



dx dy 

where, the last component in the right hand side has been drawn from the Equation [31 



(14) 



In a force- free condition (j Chandr asekhaii 1 1 96 ll : iParkerl 1 1 9 79l ) . the tension force balances 
the magnetic pressure gradient. In that condition, the Equation of magneto-hydrostatic 
equihbrium reduces to hydrostatic equilibrium 

Vp = pg (15) 

which represents a state of lowest order pressure and density. In the cases when B is strong, 
the departures from the lowest order hydrostatic pressure and density will be large. Thus 
a small variations in B can lead to large variations in plasma parameters, and could be 



one explanation for existence of the variety of fine structures ( ISu et al.l 



2009 



Tiwari et al. 



2009b 



Tiwari 



20091 ) in the strong magnetic field system of sunspots as observed from the 



high resolution data of SOT/SP. 

In the Equation [T51 the plasma pressure scale height will be independent of the scale 
height of the magnetic field. Whereas, in a tension-free field, the z-component of the 
magneto-hydrostatic equilibrium requires. 



d/dz{B^/8n +p) = pg. 



(16) 



The field in this condition cannot be force-free. In the case when tension is zero, the 
scale height of the magnetic field becomes very large, and becomes vulnerable to magnetic 
reconnection with the ov er -lying pre-existing field syste m. This can lead to flare initiation 



( 1 Venkat akr ishnan 



1990a 



Venkat akr ishnan et al. 



19931 ). 



We have computed tension force using Equation HH and expressed it in the units of 
dynes/cm^. It was found in the analysis of two active regions namely NOAA ARs 10933 
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and 10930 by 



Venkatakrishnan &: Tiwaril (120101 ) that the magnitude of vertical tension force 



attains values comparable to the force of gravity at several places over the sunspots. This 
must mean that the non-magnetic forces will not be able to balance this tension force. 
Only gradient of the magnetic pressure can match this force. These sunspot features where 
the magnetic tension balances the magnetic pressure gradient result into the force-free 
configurations. This serves as sufficient condition for verifying the force-freeness of the 
sunspot magnetic fields. We have evaluated the vertical tension force over all the sunspots 
that we have studied as listed in Table 1. 



3. Data Sets Used and Its Analysis 



We have used the high spatial resolution vector magnetograms of 19 acti ve regions (AR s) 



obtain e d from the Solar Optic a. 



(120081); 



Suematsu et a. 



Hinode (IKosugi et al. 



1120081); 



Telescope/Spectro-p o 



Ichimoto et al 



(120081 ): 



arimeter (SOT/SP: 



Shimizu et al. 



Tsuneta et al. 



( 120081 1) onboard the 



20071 1. A series of vector magnetograms of active regions (ARs) 



NOAA 10930, 10933, 10960 and 10961 have also been used to investigate the temporal 
changes in the force-free nature of the sunspot magnetic fields and also to find any 
relationship with the associated fiare activities. 

The calibration of the Hinode (SOT/SP) data have been done by using the standard 
"SP_PREP" routine developed by B. Lites which is available in the Solar Software package. 
The "SP_PREP" routine first computes the thermal shifts in the spectral and slit dimensions, 
and then applies the drift corrections for calibrating the data from levelO to levell. The 
prepared polarization spectra have then been inverted to obtain vector magnetic field 



components using an Unno-Rachkowsky (lUnno 



1956 



Rachkowskv 



the assumption of Milne- 



1982 



Skumanich &: Lites 



1967) inversion under 



ddington (ME) atmosphere (iLandolfi fc Landi Degl'Innocenti 



19871 ). The ME code used for inversion is the "STOKESFIT" 
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which have been developed by T. R. Metcalf and has been kindly made available in 
the Solar-Software package. We have used the latest version of the inversion code 
which returns the true field strengths along with the filling factor. Some of the data 
sets are also used from the Community Spectro-polarimetric Analysis Center (CSAC; 



http : / /www ■ csac . hao . ucar . edu/ ) . 



Due to insensitivity of Zeeman effect to orientation of the transverse fields, an inherent 
180° ambiguity in the azimuth determination will be present. Numerous techniques have 
been developed and applied to resolve this problem, but still complete resolution is not 
possible. The chromospheric and coronal structures are also proposed to complement the 



other methods. In the data sets which we have st udied, t 



have been removed by using acute angle method (IHarvey 



Cuperman et al. 



le 180 ° azimuthal ambigu ity 



1969 



Sakurai et al. 



1985 



19921 ). Most of the data sets used, have high spatial sampling with ~ 0.3 
arcsec/pixel. A few samples are also observed in "Normal mode" of SOT with a spatial 
sampling of ~ 0.16 arcsec/ pixel. To minimize the proiecti on effects, the vector fields are 



transformed to disk center fiVenkatakrishnan fc Gary 



19891 ) whenever they are much away 



from disk center, as is indicated in Tables 1 & 2. 

We minimize the noise present in the data sets in the following way: the pixels 
with transverse [Bt) and longitudinal magnetic field (5^) values greater than a certain 
level are only analysed. A quiet Sun region is selected for each sunspot to decide this 
critical threshold. Then la deviation in the three vector field components 5^, By and 
Bz are evaluated separately. The resultant deviation in B^ and By is taken as la noise 
level for transverse field components. Those pixels with longitudinal and transverse fields 
simultaneously greater than twice the above mentioned noise levels are only analysed. Thi s 



method to minimiz e the noise level has been successfully practiced (see , e.g. 



2009b 



Tiwari 



2009 



Tiwari et al. 



2010 



Venkatakrishnan fc Tiwari 



2009 



2010 



iwari et al. 



Gosain et al. 
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20091 . l2010l etc.). 



In most of the active regions studied, the flux is balanced with less than or equal to 
10% uncertainties. But some sunspots are single polarities with much larger imbalance of 
flux. Inclusion of such sunspots are justified as follows: we compare the values of F^/Fp, 
Fyl Fp and F^/Fp by taking whole active region and each polarity separately for those cases 
when both polarities are available. It is found that when the flux is balanced i.e., both 
polarities are taken together, the values of Fx/Fp, Py/Pp F^/Fp are smaller than the 
values computed by taking the positive and negative polarities separately. Thus, in the 
cases when flux is highly imbalanced, we consider that the lower limits of the components 
of the net Lorentz forces are obtained. It means they are more force-free than they appear 
from the calculated values. Further, the flux balance is not a concern for sufficient condition 
as we compute the vertical tension force locally. 



4. Results 



We have studied the characteristics of Lorentz forces in the photospheric magnetic fields 
of 19 sunspots obtained from the high spatial resolution spectro-polarimetric observations 
by the instrument SOT/SP aboard the Hinode satellite. The values of Lorentz force 
components normalized with the respective magnetic pressure forces are given in Table 
1. It can be noticed that the values of \Fx/Fp\, \Fy/Fp\ and \Fz/Fp\ in most of the cases 
are approximately equal to 0.1 or, sometimes even smaller. In some cases when they are 



greater than 0.1, they are always lesser than 0.6. If the ratios | F^. / F„ 



are less than unity then the field is considered to be force-free (iLow 



, \FJ FJ and \FJF 



1985 



Metcalf et al. 



1995 



Moon et al. 



2OO2I ). Thus, we can conclude that the sunspot magnetic fields are not 



far from the force-free nature as has been suspected for long time. In contrast they show 
much closeness to the approximation of the force-free atmosphere in the photosphere. Our 
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results are consistent with that of iMoon et al.l (120021 ). however, the quahty of the data used 
in our analysis is far better. 

We can notice that F^/Fp is negative for all the sunspots, except for one. The reason 
for negative Fz/ Fp can be understood in the following way: The plasma pressure is weaker 
in magnetic field concentrations such as in the sunspots as compared to the quiet Sun. 
Thus the gravity alone cannot cope with the upward pressure gradient force over sunspots. 
Therefore, to maintain the magnetohydrostatic equilibrium, downward Lorentz force [F^) 
is needed to balance excessive upward pressure gradient force. However F^ can be positive 
or zero if the flux tubes lie very low in the atmosphere. There is yet another plausible 
explanation for the negative F^/Fp. The sunspot magnetic field is confined by strong 
exterior photospheric pressure. Since the sunspot field expands outward with the height, 
the Wilson effect would be equivalent to the expanded field pushing the photospheric 
plasma downward around the edge of the sunspot. If this effect dominates, then F^/Fp will 
be negative in the layer of the atmosphere where the polarimetric signals originate. On the 
other hand, if the field is loaded with the mass in the penumbral region, its weight can push 
the field downward giving rise to a positive F^/Fp to support that weight. If this depressed 
part of the field is within the layer of polarimetric measurement, we observe positive F^jFp 
(B. C. Low 2011, private communication). 

As discussed earlier, the above results are drawn merely from the necessary conditions 
as these reflect only global property of the sunspot magnetic fields. However, as sufficient 
condition for the force-freeness, we need to look at the distribution of Lorentz forces in the 
fine structures of sunspots. We computed the vertical component of the magnetic tension 
force, a component of the Lorentz force, by using Equation [TH For example, the vertical 
tension force distribution for two sunspots, one quiet and one complex, are shown in the 
Figures 1 and 2 respectively. The associated histograms are plotted in Figures 1(c) and 2(c) 
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respectively. We distinguish between the nature of the vertical tension forces within the 
umbra and the penumbra, and this has been described later in Section 4.1. 



4.1. Nature of magnetic field in umbra and penumbra 

It can be noticed that the magnitude of vertical tension at many places (almost 
everywhere in umbra and inner penumbra) reaches up to 10~^ dyn/cm^. This value of the 
tension force is comparable with the gravitational forces experienced over sunspots which 
ranges from 5 x 10^^ dyn/cm^ to 10^^ dyn/cm^ for quiet sun density to sunspot umbral 
densities. Such high values of tension force can not be balanced only by non-magnetic 
forces i.e., gas pressure force and gravitational force. For example, when one component of 
the lorentz force (say tension force) becomes comparable to gravity, and other component 
(magnetic pressure force) is not matching, then the gas pressure force needs to become 
double to balance both gravity and magnetic tension. This being unlikely, we must only 
rely on magnetic pressure gradient to balance it. For an estimation of magnetic pressure 
gradient, let us make use of typical field strength val ues of ~ 2000 G and ~ 1000 G for 



Borrero et al 



20081 ) . Considering 



spines and intraspines respectively in penumbra (e.g. 
a vertical scale height A of ~ 100 km for magnetic field variation, we can estimate the 
magnetic pressure gradient [AB"^ /SnA) which comes out to be ~ 1.2 x 10^^ dyn/cm^. The 
magnetic tension force achieved in our study is comparable to this value at most of the 
places over sunspots. Thus, at those places where tension force is high, it is balanced by 
gradient of magnetic pressure leading to a force-free configuration. However, one should 
remember that the magnetic fields do not show a force-free behaviour over whole sunspot. 

The contours shown in Figures 1(a), 1(b), 2(a) and 2(b) distinguish between the 
umbra and the penumbra of these two sunspots. The histograms of values in the umbra 
and penumbra are plotted separately in the panels 'd' of each figure. The blue and black 
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Fig. 1. — Continuum intensity (a) and vertical tension force (T^) map (b) of a typical quiet 
sunspot NOAA AR 10933, observed from Hinode in normal scan mode of SOT/SP on 05 
January 2007 between 12 to 13 UT. We can note that the magnitudes of are high at most 
of the places over sunspot. For better contrast, the magnitudes of Tj ranging between -0.01 
and 0.01 are only plotted. However the histogram of full distribution of magnitudes of 
are shown in the panel (c). The white and black colored contours represent the boundaries 
of the umbra and the penumbra respectively. Histograms of in the umbra and penumbra 
are shown separately in the panel (d). The black and blue colored contours are for values 
in the umbra and the penumbra respectively. We can infer from the peaks of the histograms 
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Fig. 2. — Continuum intensity (a) and vertical tension map (b) of a typical active/complex 
sunspot NOAA AR 10930, observed from Hinode (SOT/SP) on 11 December 2006 between 
13 to 14 UT. We can note that the magnitudes of are high at most of the places over 
sunspot. For better contrast, the magnitudes of ranging between -0.01 and 0.01 are only 
plotted. However the histogram of full distribution of magnitudes of Tz are shown in the 
panel (c). The white and black colored contours represent the boundaries of the umbra 
and the penumbra respectively. Histograms of in the umbra and penumbra are shown 
separately in the panel (d). The black and blue colored contours are for values in the 
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coloured histograms are for values in the penumbra and umbra respectively. It is clearly 
evident from the peaks of both the histograms that the magnitudes of are higher in the 
umbra in comparison to the penumbra, thus indicating that the umbral fields are more 
force-free than the penumbral fields. From images, we can further note that the inner 
penumbra also shows high magnitudes of T^, thus indicating its force- free nature. The 
result that the umbral and inner penumbral fields are more force-free in comparison to 



the midd 
(see, e.g. 



e and the outer penumbral fields is in a greement with some earlier observations 



Solanki et al 



1993 



Mathew et al 



20041 ). The overall penumbral field seems 



to bear a mixed nature of magnetic fields depending upon nature of its fine structures 
(penumbral fibrils). Stronger ten sion forces are found in stronger and more vertical fields 



( IVenkatakrishnan fc Tiwari 



2010l ). Therefore, we can also conclude that, in the penumbra. 



the fields are more force-free in the spines which are more vertical and stronger fields. 

We have examined the behaviour of the tension force in all the active regions listed in 
Table 1. Similar behaviour is found in all the sunspots. The active regions at most of the 
places (mainly in the umbra and in the inner penumbra) show force-free behaviour. This 
serves as sufficient condition for the validity of force-freeness of the sunspot magnetic fields. 
Thus the approximation of sunspot fields as force-free magnetic fields for several purposes, 
should give reliable results. 

The mean and standard errors of Tz are given in the Tables 1 and 2. The standard 
error is computed as: 

/^T. = ^ (17) 

where is the standard deviation of vertical tension force distribution and n is the total 
number of pixels used in the computation. We can note that the average of T^ is always 
negative. This is in agreement with the magnetic fiux rope models for sunspots (see, e.g. 



Chitre 



1963 



Meyer et al. 



19771 1 which mainly dealt with the problem of global equilibrium 
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of sunspots. 



4.2. Non-linear behaviour of the force-free fields 

Once known that the sunspot magnetic fields are nearly force-free, we proceed to check 
whether the linear or non-linear force-free field approximation will be closer to validity. 
For this purpose, we examined the distribution of the force-free parameter a over all the 
active regions that are listed in Table 1. Two examples of a distribution over active regions 
are shown in the left panels of Figure 3. The NOAA AR 10933 is a fiare quiet and 10930 
is a fiare productive sunspot. Right panels show the corresponding histograms of local a 
distribution over these active regions. We find, unsurprisingly, that the sunspots with very 
complex structures such as NOAA AR 10930 show non-linear force-free field behaviour. 
Surprisingly, in the similar way, simple and quiet sunspots like NOAA AR 10933 also show 
non-linear force-free field behaviour. However the standard deviations are smaller in such 
simple active regions as compared to the complex active regions. This can be noticed, as 
an example, from the histograms in upper and lower right panels of Figure 3. We can 
observe that the range of standard deviations in both the cases are different, however, the 
distributions show similar pattern. This leads us to believe on the validity of non-linear 
force-free field approximation for both the simple and the complex sunspots. 



There are several methods proposed for computation of global 



Pevtsov et al. 



2009a 



Tiwari 



1995 



Leka fc Skumanich 



1999 



Hagino fc Sakurai 



alpha values (see, 



2004J : 



Tiwari et al. 



20091 . etc.). Most of these methods weigh on magnetic field s trength with 



the aim to predict the fiare activity by computing global twist of sunspots (iTiwari et al 



2009al ). However in the present analysis, we are merely interested in the nature of linear or 
non-linear distribution of a. Thus we compute a locally over an active region and take its 
simple mean as a global value of a for that AR. The mean and standard errors of a for all 
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the vector magnetograms are given in both the Tables. The standard error in computation 
of a has been derived from: 

fla = — (1^ 



n 



where is the standard deviation of a distribution and n is the total number of pixels 
used in the computation. More detailed statistical investigations of distribution of a can be 



found elsewhere (see, e. g. 



Pevtsov et al. 



Hagino fc Sakurai 



1994; 



Leka fc Skumanich 



1999 



Bleybel et al. 



2002 



2004 etc.). 



4.3. Temporal evolution of sunspot fields 

Further, we studied temporal variation of the force-free nature of the sunspot 
magnetic fields. For this purpose, we have selected 4 sunspots for which total of 60 vector 
magnetograms are taken from SOT/SP. Table 2 gives details of the active regions analysed 
for this purpose. The components of the net Lorentz force are plotted in Figure 4. The 
results show that there are variations in the three components with time but they are 
always less than unity. Most of the time, the ratios F^, Fy and Fz are close to 0.1 and 
less than 0.4. However in a few cases the magnitudes go beyond at certain stages of their 
temporal evolution. In those cases also the magnitudes of F^., Fy and F^ are always less 
than 0.6. Thus, the fields never obtain a magnetic field configuration that is much away 
from the force-free approximation. The distribution of vertical tension force also shows 
consistency and doesn't vary drastically with time. 

From Figure 4, we note that F^/Fp, Py/Pp F^/Fp are mostly close to zero in 
the cases of NOAA ARs 10930, 10933 and 10961 and never go beyond the value of 0.4. 
However, it is noticeable that F^/Fp, FylFp and F^/Fp goes up to around 0.5 at certain 
stages during their temporal evolution for NOAA AR 10960. This might be indicative 
of some relationship between the flare activity and the force-free field nature of sunspots. 
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a Map of NOAA AR 1 0933(05Jan07:1 21 3UT) 



400 



300 



200 



100 



: 1 ^ ' 1 1 ' 


-.■ j^f' ■ -: 













100 

E 

CO CO 

o I 

8 

-50 
-100 



100 



200 300 
pixels 



400 



a Map of NOAA AR 1 0930(1 1 Dec06:1 31 4UT) 
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Fig. 3. — Maps of the force-free parameter a for the two active regions 10933 and 10930 
respectively. We can note that the a is distributed non-hnearly over both the sunspots. 
Histograms of the a distribution of the NOAA ARs 10933 (quiet) and 10930 (flare productive) 
are shown in right panels respectively. We can see that the a is not focussed on a very close 
value in either case. Although the standard deviation is approximately three times higher in 
the case of highly flare productive sunspot 10930 in comparison to the simple sunspot 10933, 
both the sunspots show a tendency to evince non-linear force-free field atmosphere. 
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At this stage, however, we are unable to decide on this with the hmited number of the 
flare-productive sunspots and in want of their time series vector magnetograms. 

The vertical lines in the plots of Figure 4 represent the timings of the associated 
flares. The dashed line represents the timing of X-class flare (only one on 13 Dec 2006), 
dashed-dotted lines represent timings of the C-class flares and dotted lines represent timings 
of the B-class flares. As can be noted from the plots that, even if Fx/Fp, Py/Pp ^i-nd F^/Fp 
depend on the evolutionary stages of the sunspot magnetic flelds, we do not flnd any 
systematic relationship between the force-freeness of the sunspot flelds and the associated 
flare activity. 

The computation of fleld gradients might be influenced by the optical depth corrugation 
effects. An estimation however shows that the relative variation in the fleld strength and 
thus in the tension force will be 1 0%, which is much smaller thaii the observed spatial 



variation of vertical tension force (jVenkatakrishnan &: Tiwari 



20101). Acute angle method 



for azimuthal ambiguity resolution works well whenever transverse magnetic fleld vector 
makes an angle of less than 90 degrees with potential transverse fleld. Most of the active 
regions studied are simple sunspots with no complex polarity inversion lines except in the 
NOAA AR 10930. Thus, ambiguity is resolved properly in most of the active regions. We 
are interested in the magnitude of vertical tension force in the present study and not in the 
sign of it, thus we can ignore the influence of ambiguity disambiguation at neutral lines 
whenever present. Another possibility of error is because of the large scanning time of 
Hinode (SOT/SP). If the evolution of magnetic fleld is faster than the scan time of SP then 
this might cause an error in the results. This is unlikely because of larger evolution time of 
sunspot magnetic flelds, unless the sunspot is involved in an ongoing flare. For an example 
if we consider 60 arcsec size of an active region and Alfven speed of about 10 km/sec on the 
photosphere, then the Alfven travel time across the active region will be of the order of 1 
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Fig. 4.— Plots of temporal evolution of F^/Fp, Fy/Fp and F,/Fp for NOAA ARs 10930, 
10933, 10960 and 10961 as given in Table 2. We can note that the magnitudes are always less 
than unity, mostly close to 0.1 and never go beyond 0.4 except for AR 10960, where it goes 
up to 0.55. The vertical dotted lines represent timings of B-class flares, dashed-dotted hues 
represent timings of C-class flares and dashed hne (only on Dec 13 2006) represents X-class 
flare. Wc can see that, although F^/Fp, Fy/Fp and F^/Fp depend on the time evolution. 
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hour which is not smaller than the scan time of SP for such field of view. 



Discussion and Conclusions 



Using high spatial resolution spectro-polarimetric data obtained from SOT/SP aboard 
the Hinode, we find that the photospheric sunspot magnetic fields can be approximated 
to the force-free field configurations. The necessary and the sufficient, both conditions 
dealing with global and local nature of sunspot magnetic fields respectively, indicate 



that the sunspot m agnetic field s are n early force- 



results obtained by 



Moon et al. 



ree. This result is consistent with the 



Moon et al. 



(120021 ) verified the force-free 



(]2002f ). However, 

approximation by using only necessary conditions given by inequalities in Equation 10. 
We have also investigated the spatial distribution of vertical tension force which serves as 
sufficient condition. After examining the distribution of the magnitude of vertical tension 
force over all the sunspots listed in Table 1, we conclude that the magnetic fields in 
sunspot fine structures over most of the sunspot area are nearly force-free. We, further, 
distinguish between the nature of the umbral and the penumbral magnetic fields. We find 
that the umbral fields are more force-free than the penumbral fields. We also find that over 
penumbra, inner penumbra are more force free than the middle and outer penumbra. From 
the distribution of vertical tension force (Tj), it is found that t he stronger and more vertica l 



fields (spines) exhibit higher magnitudes of (cf. Figure 2 of IVenkatakrishnan fc Tiwari 



(120101 )). Thus, we can also conclude that the spines are more force-free in the penumbra. 

I 



Firstly, a study with similar purpose was performed by iMetcalf et al. 



(119951 ) using 



observations of the Stokes parameters at six wavelengths within the spectral lines of 
Na I A5858 A using derivative method. They inspected the nature of Lorentz forces by 
using the necessary conditions. Their results that the magnetic field is not force-free in 
the photosphere and becomes force-free above certain heights in the lower chromosphere 
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has been subje cted to an undere stimation of transverse field due to their weak-field 



approximation ( iMoon et al. 



2OO2I) 



The calculations of vertical current density and alpha param e ter over sunspots obs erved 



from SOT/S 



Su et al. 



mm; 



Tiwari et al. 



fl2009bf l: 



P has been performed recently by 
Tiwaril ( l2009l ). The distribution of a is found to be non- uniform over sunspots. After 
verifying that the sunspot fields are almost force-free, we have examined the range of the 
distribution of alpha over several sunspots, including simple and complex active regions, to 
inspect the imminence of linear or non-linear force-free field approximation. Our results 
show that the non-linear force-free field approximation is applicable not only in complex 
active regions such as NOAA AR 10930 but also in the simple active re gions such as NOA A 



AR 10933 (see Figure 3). This result is not in accordance with that of 



Moon et al. 



(12002[ l. 



Th e reason cou l d be t he difference in the quality of data used. The resolution of data used 



by 



Moon et al. 



(|2002[ ) might have not been sufficient to obtain the non-linear behaviour 
of alpha over sunspots they studied. Thus, the investigations using high resolution data 
lead us to believe that the sunspot magnetic fields are nearly force-free in nature and they 
imitate non-linear force-free behaviour. 

A study of temporal evolution of Lorentz force components in 60 vector magnetograms 
of four active regions (see Table 2 and Figure 4) shows variations in the Lorentz force 
components normalized by Fp i.e., \Fx/Fp\, \Fy/Fp\ and \Fz/Fp\. They are always lesser than 
unity, mostly closer to 0.1 and never show higher values greater than 0.6, thus remaining 
under the validity of the force-free field approximation. These sunspots, in every phase of 
their evolution, show immediacy with the force-free approximation of their magnetic fields. 
From Figure 4, we further tried to find any relationship between the temporal evolution of 
the force-free field nature of sunspots and the associated fiare- activity. Although \Fx/Fp\, 
\Fy/Fp\ and \Fz/Fp\ depend on the evolutionary stages of active regions, they don't show a 
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systematic relationship with the occurrence of flares associated with those active regions. 



We depend on the extrapolations of photospheric magnetic fields to obtain the coronal 
vector magnetic fields owing to inadequacy of coronal field measurements. As it is well 
known and also discussed in the introduction, that all the extrapolation techniques have to 
rely on the photospheric vector field measurements and also on the force-free approximation 
of these magnetic fields. The effect of several possible inaccuracies in the vector field 
measurement on the n on-linear force-free field modelings has recently been studied by 



De Rosa et al. 



mm; 



Wiegelmann et al. 



feOlOl ). Coronal magnetic field reconstruction 



by extrapolations of photospheric magnetic fields under non-linear m odelings have shown 



satis: 



' actory results by nearly matching with the coronal observations (IMcClvmont fc Mikic 



1994 



Bleybel et al. 



2002 



Wiegelmann et al. 



2005 



Schrijver et al. 



2006 



Met calf et al. 



20081 ). These results then also support our conclusion that the sunspot magnetic fields are 



close to satisfy the non-linear force-free field approximation. 



A preprocessing procedure has been developed by 



Wiegelmann et al. 



(jaOOd) to make 



the observed data completely force-free as a suitable boundary conditions for force-free 



extrapo^ 
used by 



ations. This has led mor e reliable better extrapolations of the fields. The method 



Wiegelmann et al. 



( 120061 ) however uses only necessary condition to make the fields 
force-free. Also, most of the time observed field of view suffers from the flux imbalance 
leading to an inaccuracy in the preprocessing. The preprocessing might be improved by 
using the sufficient condition (as described in Section 2.2) in preparing the boundary data 
for the force-free fleld extrapolations. 

There is, however, good enough scope of improvement of the present study. At present 
our investigation is limited to comp utation of vertical c ompo nent of the magnetic tension 



force. Whereas in a recent study by 



Puschmann et al. 



(l2010l ). it has been found that the 



horizontal component of Lorentz force is greater than the vertical component. Detailed 
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analysis of full sunspot magnetic fields within certain optical heights are required to 
understand full scenario of force-free nature of the sunspot magnetic fields. In a forthcoming 
work, we plan to look into all the three components of the equilibrium forces in the height 
dependent inversions of a complete sunspot. 
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1. List of the Active Regions, that we have studied, with their observational details and resiilts obtained. 



NOAA AR 
Number 


Date: Time(UT) 
of Observation 


Position 


F3:/ Fp 


Fy/Fp 


FJFp 


< O > ±Ma 

Xl0-*m-l 


< 
XlO 


> ±MT^ 


comment 


10972 


07 Oct 2007: 


0200 


S05W10 


-0.137 


0.203 


-0. 


.600 


0.225 ± 0.819 


-6, 


.529 ± 0, 


.067 


quiet, diffused 


10971 


29 Sep 2007: 


1011 


N03W04 


0.006 


0.102 


-0. 


.608 


0.648 ± 0.768 


-6, 


.038 ± 0, 


.063 


quiet, diffused 


10970 


05 Sep 2007: 


0203 


S07W46(t) 


-0.155 


0.213 


-0 


.148 


2.655 ± 0.731 


-3 


.446 ± 


,047 


quiet, diffused 


10969 


29 Aug 2007: 


0001 


S05W33(t) 


0.384 


0.131 


-0, 


.253 


2.930 ± 0.698 


-5, 


.735 ± 0, 


,049 


quiet 


10966 


07 Aug 2007: 


0500 


S06E20(t) 


-0.028 


0.196 


-0, 


.346 


-2.249 ± 1.364 


-6, 


.303 ± 0, 


,086 


complex, weak 


10963 


12 Jul 2007: 


0001 


S06E20(t) 


-0.003 


0.119 


-0, 


.206 


-3.252 ± 0.377 


-4, 


.730 ± 0, 


,031 


active 


10961 


01 Jul 2007: 


0304 


S13E06(t) 


0.024 


0.172 


-0, 


.271 


-4.754 ± 0.939 


-6, 


.753 ± 0, 


,072 


quiet 


10960 


07 Jun 2007: 


0304 


S07E07 


0.137 


0.093 


-0, 


.482 


-1.588 ± 0.446 


-6, 


.625 ± 0, 


,046 


complex, active 


10956 


18 May 2007: 


0809 


N03W01 


0.248 


-0.037 


-0, 


.237 


8.523 ± 0.636 


-6, 


.974 ± 0, 


,060 


complex, rotating 


10955 


13 May 2007: 


1100 


S09W28(t) 


0.254 


-0.035 


-0, 


.349 


0.411 ± 0.689 


-6, 


.143 ± 0, 


,063 


quiet 


10953 


29 Apr 2007: 


0001 


S10E22(t) 


0.046 


0.067 


-0, 


,089 


-2.499 ± 0.609 


-4 


.424 ± 


,058 


quiet, strong 


10944 


03 Mar 2007: 


0001 


S05W30(t) 


0.078 


-0.032 


-0, 


.133 


-1.749 ± 0.284 


-3, 


.731 ± 0, 


,029 


quiet, strong 


10940 


01 Feb 2007: 


0708 


S04W06 


0.055 


-0.033 


-0, 


.279 


-1.265 ± 0.279 


-7, 


.371 ± 0, 


,032 


complex 


10939 


22 Jan 2007: 


1314 


S04W46(t) 


0.148 


-0.034 


-0, 


.086 


-1.317 ± 0.660 


-6, 


.512 ± 0, 


,067 


complex 


10933 


05 Jan 2007: 


1213 


S04E03 


0.045 


0.002 


-0, 


.108 


-0.557 ± 0.123 


-6, 


.341 ± 0, 


,026 


quiet 


10930 


11 Doc 2006: 


1314 


S05W()7 


0.039 


0.038 


-0, 


,062 


-8.604 ± 0.316 




.749 ± 


,028 


complex, rotating 


10926 


03 Dec 2006: 


0607 


S09W26(t) 


0.139 


0.114 


-0, 


.183 


1.674 ± 0.849 


-5, 


.479 ± 0, 


,068 


complex 


10923 


10 Nov 2006: 


1617 


S06E40(t) 


0.152 


0.007 


0.037 


-4.840 ± 0.761 


-4, 


.903 ± 0, 


,069 


complex, strong 


10921 


06 Nov 2006: 


1416 


S08W36(t) 


0.195 


0.121 


-0.222 


-1.169 ± 0.463 


-4, 


.603 ± 0, 


,041 


complex 


(t) : transformed 
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Table 2: noaa 

obtained. Associated 



ARs 10930, 10933, 10960 and 10961 
plots of F^/Fp, Fy/Fp and F^/Fp 



during tiicir temporal evolution for few days with their observational details and results 
are shown in Figure 4. 





Date: 


Time(UT) 


Position 


Fx/ F'p 


Fy/Fp 


Fz/Fp 


< a > ±^c, 


<Tz > ±MT, 


comment 


Number 


of Observation 












XlO~^m" 


-1 


XlO 


^dyn 


cm ^ 






10930 


10 


Dec 


2006: 


0102 


S04E14(t) 


0.078 


-0.003 


-0. 


,069 


-2.854 ±0.719 


-6. 


,601 


± 


0, 


,061 


complex, 


rotating 


10930 


10 


Dec 


2006: 


2122 


S04E03 


0.115 


0.022 


-0. 


,051 


-6.691 ± 0.699 


-7, 


,609 


± 


0, 


,048 


complex, 


rotating 


10930 


11 


Dec 


2006: 


0304 


S04W01 


0.079 


0.029 


-0, 


,076 


-7.875 ± 0.536 


-7, 


,267 


± 


0, 


,045 


complex, 


rotating 


10930 


11 


Dec 


2006: 


0809 


S04W04 


0.059 


0.033 


-0 


089 


-6.946 ± 


,508 


-6, 


,213 


± 


0, 


044 


complex, 


rotating 


10930 


11 


Dec 


2006: 


1112 


S04W06 


0.047 


0.034 


-0, 


,086 


-10.252 ± 0.534 


-6, 


,462 


± 


0, 


,048 


complex. 


rotating 


10930 


11 


Dec 


2006: 


1314 


S05W07 


0.039 


0.038 


-0. 


,062 


-8.604 ± 0.316 


-5, 


,749 


± 


0, 


,028 


complex. 


rotating 


10930 


11 


Dec 


2006: 


1718 


S05W08 


0.016 


0.044 


-0. 


,077 


-9.617 ± 0.662 


-6, 


,636 


± 


0, 


,050 


complex. 


rotating 


10930 


11 


Dec 


2006: 


2021 


S05W10 


-0.001 


0.047 


-0, 


,088 


-10.035 ± 0.555 


-6, 


,617 


± 


0, 


,051 


complex. 


rotating 


10930 


11 


Dec 


2006: 


2324 


S05W12 


-0.013 


0.045 


-0 


063 


-10.490 ± (1 


1.595 


-6 


,401 


± 


0, 


059 


complex. 


rotating 


10930 


12 


Dec 


2006: 


0304 


S06W14(t) 


0.099 


0.039 


-0. 


,029 


-13.988 ± 0.599 


-5, 


,959 


± 


0, 


,069 


complex. 


rotating 


10930 


12 


Dec 


2006: 


1011 


S06W18(t) 


0.102 


0.062 


-0. 


,024 


-11.561 ± 0.325 


-4, 


,971 


± 


0, 


,047 


complex. 


rotating 


10930 


12 


Dec 


2006: 


1516 


S06W21(t) 


0.119 


-0.060 


-0, 


042 


-8.919 ± 


,547 


-5, 


,261 


± 


0, 


,056 


complex. 


rotating 


10930 


12 


Dec 


2006: 


1718 


S06W22(t) 


0.118 


-0.063 


-0, 


,044 


-8.342 ± 0.509 


-5, 


,359 


± 


0, 


,053 


complex. 


rotating 


10930 


12 


Dec 


2006: 


2021 


S06W24(t) 


0.127 


0.066 


-0. 


,047 


-8.803 ± 0.639 


-5, 


,697 


± 


0, 


,063 


complex. 


rotating 


10930 


13 


Dec 


2006: 


0406 


S06W29(t) 


0.137 


0.078 


-0. 


,025 


-11.117 ± 0.551 


-6, 


,389 


± 


0, 


,069 


complex, 


rotating 


10930 


13 


Dec 


2006: 


0708 


S06W31(t) 


0.136 


0.079 


-0, 


,002 


-9.793 ± 0, 


,617 


-6, 


,153 


± 


0, 


,076 


complex. 


rotating 


10930 


13 


Dec 


2006: 


1314 


S06W34(t) 


0.125 


0.083 


-0, 


,007 


-6.289 ± 0, 


,302 


-5, 


,001 


± 


0, 


,034 


complex. 


rotating 


10930 


13 


Dec 


2006: 


1617 


S06W35(t) 


0.126 


0.084 


-0. 


,006 


-7.712 ± 0, 


,555 


-6, 


,061 


± 


0, 


,060 


complex. 


rotating 


10933 


04 


Jan 


2007: 


1112 


S04E19(t) 


0.202 


0.029 


-0. 


,102 


-1.082 ± 0, 


,181 


-6, 


,266 


± 


0, 


,026 


quiet 


10933 


04 


Jan 


2007: 


1819 


S04E15(t) 


-0.049 


0.018 


-0 


174 


-3.019 ± 0, 


,308 


-6, 


,110 


± 


0, 


,047 


qui 


et 


10933 


05 


Jan 


2007: 


0001 


S04E11 


0.125 


0.017 


-0, 


,146 


-1.254 ± 0, 


,233 


-6, 


,721 


± 


0, 


,050 


quiet 


10933 


05 


Jan 


2007: 


0708 


S04E06 


0.086 


0.003 


-0. 


,113 


-0.569 ± 0, 


,233 


-6. 


,624 


± 


0, 


,043 


quiet 


10933 


05 


Jan 


2007: 


1213 


S04E03 


0.046 


0.002 


-0. 


,108 


-0.557 ± 0, 


,123 


-6, 


,341 


± 


0, 


,026 


quiet 


10933 


06 


Jan 


2007: 


0001 


S03W03 


-0.032 


-0.017 


-0, 


,134 


-0.520 ± 0, 


,170 


-5, 


,765 


± 


0, 


,038 


quiet 


10933 


06 


Jan 


2007: 


0405 


S03W06 


-0.065 


-0.018 


-0, 


128 


-0.406 ± 


,120 


— .5 


,650 


± 


0, 


025 


quiet 


10933 


06 


Jan 


2007: 


1000 


S03W09 


-0.106 


0.001 


-0. 


,112 


-0.847 ± 0, 


,112 


-6, 


,247 


± 


0, 


,028 


quiet 


10933 


06 


Jan 


2007: 


1213 


S03W10 


-0.133 


0.002 


-0. 


,095 


-0.458 ± 0, 


,132 


-6, 


,467 


± 


0, 


,030 


quiet 


10933 


07 


Jan 


2007: 


0001 


S04W17(t) 


0.135 


-0.005 


-0. 


,158 


-2.223 ± 0, 


,217 


-5, 


,175 


± 


0, 


,036 


qui 


et 


10933 


07 


Jan 


2007: 


1213 


S04W20(t) 


0.103 


-0.003 


-0. 


,131 


-2.064 ± 0, 


,261 


-4, 


,524 


± 


0, 


,036 


quiet 


10933 


07 


Jan 


2007: 


1819 


S04W25(t) 


0.161 


-0.018 


-0. 


,072 


-2.382 ± 0. 


.267 


-5. 


,108 


± 


0, 


,037 


quiet 


10960 


06 


Jun 


2007: 


0607 


S08E29(t) 


0.079 


0.066 


-0. 


,144 


-3.270 ± 0, 


,610 


-4, 


,188 


± 


0, 


,049 


complex 


, active 


10960 


05 


Jun 


2007: 


1112 


S08E26(t) 


0.087 


0.079 


-0. 


,222 


-2.755 ± 0, 


,439 


-3, 


,773 


± 


0, 


,037 


complex 


, active 


10960 


06 


Jun 


2007: 


0203 


S07E19(t) 


0.119 


0.074 


-0, 


,335 


-3.423 ± 0, 


,480 


-4, 


,465 


± 


0, 


,046 


complex, active 


10960 


06 


Jun 


2007: 


0708 


S07E17(t) 


0.126 


0.076 


-0. 


,375 


-2.711 ± 0, 


,490 


-4, 


,525 


± 


0, 


,048 


complex 


, active 


10960 


06 


Jun 


2007: 


1213 


S07E14(t) 


0.111 


0.089 


-0. 


,441 


-2.226 ± 0, 


,374 


-4, 


,394 


± 


0, 


,038 


complex 


, active 


10960 


06 


Jun 


2007: 


1920 


S07E10 


0.192 


0.086 


-0. 


,464 


-1.252 ± 0. 


,373 


-4, 


,867 


± 


0, 


,037 


complex 


, active 


10960 


06 


Jun 


2007: 


2223 


S07E09 


0.198 


0.092 


-0, 


,451 


-1.105 ± 0, 


,478 


-5, 


,678 


± 


0, 


,048 


complex 


, active 


10960 


07 


Jun 


2007: 


0304 


S07E07 


0.137 


0.093 


-0 


482 


-1.588 ± 


,446 


— 5 


,625 


± 


0, 


,045 


complex, active 


10960 


07 


Jun 


2007: 


1920 


S07W03 


0.008 


0.109 


-0, 


,535 


-2.302 ± 0, 


,501 


-6, 


,127 


± 


0, 


,050 


complex 


, active 


10960 


08 


Jun 


2007: 


0708 


S07W09 


-0.178 


0.100 


-0. 


,514 


-1.631 ± 0, 


,666 


-6, 


,821 


± 


0, 


,055 


complex 


, active 


10960 


08 


Jun 


2007: 


0809 


S07W10 


-0.204 


0.103 


-0, 


505 


-1.055 ± 0, 


,570 


-6 


,928 


± 


0, 


,055 


complex 


, active 


10960 


08 


Jun 


2007: 


1314 


S07W12(t) 


-0.223 


0.129 


-0, 


,553 


0.302 ± 0.418 


-6, 


,101 


± 


0, 


,039 


complex 


, active 


10960 


08 


Jun 


2007: 


1516 


S07W14(t) 


0.226 


0.086 


-0. 


,559 


-2.764 ± 0. 


,360 


-5, 


,898 


± 


0, 


,034 


complex 


, active 


10960 


08 


Jun 


2007: 


1819 


S07W16(t) 


-0.238 


0.079 


-0. 


,533 


-4.118 ± 0, 


,396 


-6. 


,120 


± 


0, 


,036 


complex 


, active 


10960 


08 


Jun 


2007: 


2021 


S07W17(t) 


-0.208 


0.073 


-0, 


,519 


-4.120 ± 0, 


,423 


-6, 


,233 


± 


0, 


042 


complex 


, active 


10961 


29 


Jun 


2007: 


0100 


SllE32(t) 


-0.053 


0.147 


-0, 


092 


-4.429 ± 


,794 


-3, 


,592 


± 


0, 


054 


quiet 


10961 


29 


Jun 


2007: 


0700 


SllE28(t) 


-0.057 


0.147 


-0, 


,126 


-5.473 ± 0, 


,819 


-4, 


,023 


± 


0, 


,055 


quiet 


10961 


29 


Jun 


2007: 


1819 


SllE21(t) 


-0.026 


0.144 


-0. 


,176 


-4.908 ± 0, 


,876 


-4, 


,384 


± 


0, 


,061 


quiet 


10961 


30 


Jun 


2007: 


0000 


S12E19(t) 


-0.032 


0.137 


-0 


182 


-4.654 ± 


,972 


-4 


,569 


± 


0, 


073 


qui 


et 


10961 


30 


Jun 


2007: 


0405 


S12E17(t) 


-0.044 


0.149 


-0, 


217 


-6.521 ± 


,863 


-4, 


,661 


± 


0, 


062 


qui 


et 


10961 


30 


Jun 


2007: 


0910 


S12E14(t) 


-0.025 


0.167 


-0. 


,249 


-6.369 ± 0. 


,973 


-5, 


,676 


± 


0, 


,078 


quiet 


10961 


30 


Jun 


2007: 


1100 


S13E13(t) 


-0.002 


0.167 


-0. 


,307 


-4.725 ± 0, 


,616 


-4, 


,661 


± 


0, 


,039 


quiet 


10961 


30 


Jun 


2007: 


2223 


S13E08(t) 


-0.010 


0.128 


-0. 


,174 


-4.235 ± 1, 


,549 


-5, 


,579 


± 


0, 


116 


quiet 


10961 


01 


Jul 2007: 


0304 


S13E()5(t) 


0.024 


0.172 


-0, 


271 


-4.754 ± 


,939 


— 5 


,753 


± 


0, 


072 


qui 


et 


10961 


01 


Jul 2007: 


0809 


S13E02(t) 


-0.047 


0.183 


-0, 


,257 


-5.437 ± 1, 


,064 


-6, 


,404 


± 


0, 


,082 


qui 


et 


10961 


01 


Jul 2007: 


1314 


S13W01(t) 


-0.054 


0.186 


-0. 


,269 


-2.786 ± 0, 


,964 


-5, 


,879 


± 


0, 


,073 


quiet 
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